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Abstract 

The notion of operator resonances was introduced earlier by Al. Zamolodchikov within the frame- 
work of the conformal perturbation theory. The resonances are related to logarithmic divergences of 
. integrals in the perturbation expansion, and manifest themselves in poles of the correlation functions 

and form factors of local operators considered as functions of conformal dimensions. The residues of 
\ the poles can be computed by means of some operator identities. Here we study the resonances in the 

fSJ ■ Liouville, sinh- and sine-Gordon models, considered as perturbations of a massless free boson. We 

show that the well-known higher equations of motion discovered by Al. Zamolodchikov in the Liou- 
^ , ville field theory are nothing but resonance identities for some descendant operators. The resonance 

expansion in the vicinity of a resonance point provides a regularized version of the corresponding op- 
04 ' erators. We try to construct the corresponding resonance identities and resonance expansions in the 

sinh- and sine-Gordon theories. In some cases it can be done explicitly, but in most cases we are only 
able to obtain a general form so far. We show nevertheless that the resonances are perturbatively 
exact, which means that each of them only appears in a single term of the perturbation theory. 

D . 1. Introduction 

The notion of resonances in the conformal perturbation theory was introduced by Al. Zamolodchikov 
fT^ ■ in [I]- He considered a conformal field theory perturbed by a relevant perturbation, i. e. a perturbation 

1 described by the action 5pcrt ~ / (P' x ^ pcrt{x) , vifhere $port(a;) is a local spinless field of conformal 

■ dimension Aport < 1- Let 0{x) be any field of given conformal dimension A. Suppose that the fusion 
rule of the O operator and k instances of the $pert operators contains an operator O' (x) of the dimension 

■ A'. Then, if A > A'"*' = A' -I- fc(l — Apert), ultraviolet divergences take place in the fcth order of the 
perturbation theory. To avoid the divergence, the operator O must be renormalized by adding the 
operator O' with an infinite coefficient. Alternatively, the renormalization can be often achieved by 
analytic continuation from the region of convergence. The most subtle situation takes place on the 
threshold of the ultraviolet divergences A — . In this case the divergence is logarithmic and the 
definition of the renormalized field O'"^'^ is ambiguous. Another manifestation of this phenomenon is that 
every correlation function of the operator O possesses a pole at the point A ~ . The residue of the 
pole coincides with a correlation function of the operator C so that one can write the identity as an 

^ ■ operator one: 

^Res_^^ 0{x) = const xO'{x). (1.1) 

This phenomenon plays an important role in the two-dimensional integrable quantum field theory. Exact 
vacuum expectation values and form factors are given by analytic functions in the conformal dimensions 
of the corresponding operators. These functions possess poles at the resonance points 013]. In order to 
extract some physically significant information at these points, the singularities must be cancelled by an 
appropriate renormalization procedure. 

Another famous result by Al. Zamolodchikov is the discovery of the higher equations of motion in 
the Liouville field theory [3]. The Liouville theory is a theory of a scalar real boson field f{x) with the 
potential of the form e'"''. Being a conformal field theory it contains two Virasoro algebras generated by 

the modes Ck and >C/c of the two nonzero components of the energy momentum tensor T and Tj^^. 

The higher equations of motion are operator relations of the form 



X 
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Here Umn = -'-j^6~^ + -^-^b, Bmn are some constants, and D„in and Dmn are polynomials in the 
modes C-k and C-k {k > 0) correspondingly of the level mn, such that r'm„e''""''^ = Z)m„e'''""'^ = 0. 
The simplest case m ~ n = I reduces to the usual equation of motion. In the present paper the 
higher equations of motion are interpreted as resonance identities of the form (|l.ip in the perturbation 
theory from the free massless boson with the perturbation J d^x e*"^, and the corresponding rcnormalized 
operators are defined. 

The sinh-Gordon theory can be considered either as a relevant perturbation of the Liouville theory 
or as a relevant perturbation of the free massless boson. One may expect that it inherits the resonances 
of the Liouville theory, and, besides, possesses its own resonances. We shall show how the resonance 
identities are modified in the sinh-Gordon theory. It turns out that if either m or n is even the resulting 
resonance identity coincides with that in the Liouville theory, while in the case of odd m and n additional 
terms appear. These additional terms can be easily written explicitly, if either to = 1 or n = 1. Further, 
we try to generalize the results to the sine-Gordon theory. The situation in that case becomes more 
involved due to the fact that the very operator O' can possess a resonance pole. We use the form factor 
analyticity conjecture to understand the structure in this case. 

Note that one of the main motivations to this study was caused by the present unsatisfactory state 
of the form factor theory. Within the bootstrap approach we try to identify solutions to the form factor 
equations [S] with some particular operators in the Lagrangian or conformal perturbation theory, and 
ignore the fact that we do not clearly understand the very sense of these operators beyond the free field 
and conformal field theory. 

In Sec. [2] we recall some basic facts concerning the Liouville field theory and the higher equations of 
motion. In Sec. [3] we demonstrate the connection of the higher equations of motion with the resonance 
identities. The generalizations to the sinh- and sine-Gordon models are made in Sees. |3] and [5] corre- 
spondingly. In Appendix 12] the logarithmic behavior of some particular fields in the Liouville theory is 
substantiated. Appendix iBl contains some support of the results of the paper by the form factor theory. 



2. Liouville field theory and higher equations of motion 

Here we recall some basic facts concerning the Liouville field theory and, in particular, what is called the 
higher equations of motion. The Liouville field theory on the Minkowski plane is defined by the action 

5lm^ jd^x^^^^c^^y (2.1) 

We shall also use the letters Q and p defined as 



Q = b + b-\ P^Q-\ b^ = -^. (2.2) 

1 -p 



Besides, we shall always use the light-cone coordinates and derivatives: 

d 



X 



^ dx^' 

The space of local operators of the Liouville theory is spanned by the exponential operators 

Va{x) = e'"^^^) (2.3) 

and their descendants 

• • • dt"'ipd^^ip- ■ ■ d^^ifiVaix). (2.4) 

The pair of numbers [K, L) — (^ ki, ^ k) is called the level of a descendant. The numbers K and L 
are called its chiral (right and left correspondingly) levels. The descendants with L = are called right, 
while those with K = left. The (Lorentz) spin of the operator is equal to K — L. 

It is known that not all the operators (|2.4I) . and even not all the exponential operators (|2.3p . are 
different. The exponential operators are pairwise identified due to the reflection relation [B] 

Va{x) = RaVQ^aix). (2.5) 
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The function Ra is called the reflection function and given by 



r(2fo(g-a)-62)r(26-i(g-a)-6-2)' ^ 

where 

This identification imposes an identification of the descendant operators, which will be described a little 
later. 

Since these basic elements depend on the continuous parameter a, a complete description of the 
space of local operators must admit linear combinations involving integrating and derivatives in a. In 
particular, we may consider the operators 

and their descendants 9^ (/j • • • 9^:-" <^ d'^^^p ■ ■ ■ d''^ip V'^ {x) . 

The Liouville field theory possesses the conformal symmetry, which is preserved on the quantum 
level. It means the stress-energy tensor can be rendered traceless by adding some full derivatives in 

coordinates. Hence, in the light-cone coordinates we have T-| = T |^ = 0. For two other components 

we shall use the standard notation 



^gj = -27rT^^(x) = -\{d^^f + ^dl^. (2., 



Due to the energy-momentum conservation the component T is a function of x~ only, being x^-indepen- 
dent, while the component T is a function of x'^ . This defines two Virasoro algebras as follows. Consider 
the Wick rotation = — ir. The light-cone variables becomes complex coordinates on the Euclidean 
plane: x" ^ z ^ x^ + it, x'^ — z ~ x^ — ir. One can define the operators £„, £„ acting on local 
operators 0{x): 

CnO{x) = j^^{w~ zr+^T{w)0{x), UO(x) = ^ '^[w - zr+'fiw)Oix). (2.9) 
These operators form the Virasoro algebra 

Cn] = (m - n)Cm+n + -^m{m'^ - l)5m+n,0, 

and the same algebra for £„, while — 0. Both right and left algebras possess the same central 

charge 

c = 1 + 6Q2. 

With respect to the Virasoro algebras the exponential operators Va are primary, 
with the conformal dimensions 

Aa = a{Q - a). 

The descendants l|2.4p for generic values of a can be expressed as linear combinations of the conformal 
basis: 

C^k^■■■C-k,„C-l,■■■C-l,yaix), (2.10) 

where the level {K, L) coincides with (^ ki, J2 h)- Note that the reflection property for the descendant 
operators has a simple form in this basis. Namely, 

£_fc, •••£_Z„V;(x) =i?a/:-fei •••/:-fc„,^-h •••'C-/„yQ-a(2:). (2.11) 

The particular operators £_i and £_i always coincide with the light-cone derivatives: 
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Thus, the standard equation of motion 

d^d+(p = 2Trnbe'"f (2.12) 

can be rewritten as 

C-iC-iVo{x) = 27rfibVb{x). (2.13) 

In the quantum theory the equation of motion does not determine the behavior of the model as a whole. 
It just establishes a relation between several operators. Al. Zamolodchikov proposed [4, a generalization 
of the equation of motion in the Liouville field theory. Let us recall his main result. 
Consider the special type of exponential operators 

Vmn{x) ^ Va^^{x), a„i„ = ^— ^6"^ + ^-^^5, m,neZ. (2-14) 

For a = ajnn and m, n > 0, the conformal dimensions coincide with the Kac dimensions [7] and the 
module of the Virasoro algebra generated by £_/c, C-k {k > 0) on the operator Vmn{x) as on the highest 
weight vector is non-free on the level mn. It means that there exist nonzero operators of the form 



Drnn = ^ ^ ^fci . . .fc^ -^-fei ' ' ' -C-fc^ 

S — 1 l<fcl<--<fcs 
fclH hfcs=mn 

mn 

Dmn — ^fci...fcs-C-fci • • • £-fc^ 



(2.15) 



s — 1 l<ki<--<ks 

fclH hfcs="i 



with some real coefficients A..., such that 

DmnVmn = i'nmKin =0, m, n > 0. (2-16) 

It is convenient to normalize these operators by the condition = 1, so that 

n — rinn 1 . . . f) _ rrmi i . . . 

The higher equations of motion read 

Dmn-^mn^inn^'^^ '^■^ran^m^—ni,'^') • (2-17) 

Here the coefficients Bmn are explicitly known 

i?™„=A"s™„=A"&'+'"-'"-f#^^^^^^ n {kh-^+ib). (2.18) 

(fc. 0^(0,0) 

3. Higher equations of motion as resonance identities 

Let us recall the phenomenon of resonances between operators in the conformal perturbation theory [TH3] • 
Consider a field theory in two dimensions that can be represented as a perturbation of a conformal field 
theory by a relevant primary operator: 

S = Sa + 5pcrt = 5*0-9 j $port(a;), (3.1) 

where 5*0 is the action of a conformal field theory and 'Sport (2;) is a spin primary conformal operator 
of conformal dimension Aport < 1- The coupling constant g has the dimensionality (mass)^'-^~'^p"'^. 
Formally, any (chronologically ordered) correlation function (• • • ) can be represented in terms of the 
unperturbed correlation function (• • ■ )o as 

1,^ t VK \ \ (e'^^<^i(^i)g2(.T2)---)o 
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Since iSport is proportional to a "small" parameter g, the r. h. s. can be expanded in a pertm^bation series 
in g. Naturally, an arbitrary term of this expansion consists of terms of the form 

(fyi ■ --d^Vk (^-pertCj/i) ■ ■ •$pcrt(2/fc)Ci(a;i)C2(a;2) 

The integrals are calculated by their continuation to the Euclidean spacetime, so that we shall suppose 
such continuation everywhere, and write \x\'^ — —x^ = . The integrals may possess both ultraviolet 

and infrared divergences. We will be interested in the former. Consider the region of the integration 
over Di in the vicinity of the point where an operator is inserted, e. g. x — xi. Suppose that one of 
the operators in the correlation function 0{x) is any, not necessarily primary, spin operator of a given 
conformal dimension A. For simplicity, consider the case k — 1. Then applying the operator product 
expansions we obtain 

where C„ are some structure constants and A^ are conformal dimensions of the operators (x) that 
enter the operator product expansion. Operators of nonzero spin, which have different right and left 
conformal dimensions, do not contribute the expansion in the ultraviolet region. If the dimension A^ 
large enough, A^ > A + Apert — 1, the corresponding contribution is convergent. If AJ^ < A + Apcrt — 1, 
the contribution diverges, but the divergence can be canceled by renormalization of the operator O 
according to O — >■ O"" — O + a„gA^(^p'='^'+^~'^"~^^0^ with some constant factor a„ and A being the 
ultraviolet cutoff in the momentum space, which will be set to infinity after renormalization. In both 
cases the renormalized operator can be defined in a unique way. The situation when the resonance 
condition 

a; = A + Apert - 1 

is satisfied turns out to be most difficult. In this case the integral diverges logarithmically and the 
renormalization 

O'-^iX; x) = 0{x) + a,,g (^log 0'^{x) 

inevitably contains some arbitrary but finite parameter A of the dimensionality of mass. This parameter 
makes the definition of the operator 0''=" intrinsically ambiguous. The last admits a redefinition 

0""(A'; x) = 0""(A; x) - a„g (^log ^) 0'„{x). 

In terms of correlation functions the logarithmic divergence reveals itself as poles of correlation functions 
with respect to A, while the residues of the poles are proportional to the correlations functions, where 
the operator O is substituted by O^. 

Generally, consider an operator 0{x) of the dimension A and an operator 0'{x) of the dimension A' 
that appears in the expansion in the powers oi yi — x of the operator product 

C(a;)$pcrt(2/i)---*port(yfe)- 

We suppose that A is a continuous parameter^ so that the conformal dimension A' of the admissible 
operator is a function of A: A' = AJj,(A). Then every correlation function of 0{x) possesses a pole at 
the resonance point A = A,, which is a solution to the equation 

A,^Ar(A,), Ar(A)1^'A^(A) + fc(l- Apert). (3.2) 

The residue of the pole in a correlation function of O is proportional to the corresponding correlation 
function of the operator C, so that we may write it as an operator identity 

^Rej 0(A;x) = ag'=0;(x). (3.3) 



^Surely, there are examples where the set of admissible dimensions is discontinuous. In these cases one may need to 
use some tricks to stir conformal dimensions to see the poles (see e. g. [3]). Such theories are not the subject of present 
consideration. 



5 



There are two ways to define the renormahzed operator at the residue point. The first way is the 
'traditional' one by means of cutting off integrations and subtracting the divergent part: 



The parameter A can be taken, in principle, arbitrary. It means that the logarithmic divergence leads to 
an ambiguity: 



&r{\';x) ^ &r{\;x) ~ 2ag^ (^log y ) O'^. 
The other way uses the expansion in the vicinity of the pole: 



fcN2(A-Ar(A))0/^ N 

^^""^ ^ A-Ar(A) ^ ^ ' ^r(A)), (3.4) 

where A is again an arbitrary constant of the dimension of mass. The factor A^^^"^*: ^ in the first term 
is introduced to equate the dimension of the r. h. s. to that of the 1. h. s. The arbitrariness of A results 
in the ambiguity of the definition of the renormahzed operator. 

Both approaches lead to essentially the same result. Namely, there exists a constant k such that 

OT{X;x)^dr{K\-x). 

The value of A is at our disposal. If the perturbed theory is non-conformal, it possesses a mass scale 
M ~ c,i/2(i-Apert)^ and we have to take A cx M with an arbitrary dimensionless proportionality constant. 
In the conformal case we are not bound even by this limitation. 

Let us stress that the expansion (|3.4I) is taken over the powers of A — A^°^(A) rather than of A — A* , 
and the numerator of the first term contains the operator C(x), which is a function of the parameter 
A, rather than the resonance operator 0'.j,{x). This form of the expansion is dictated by the structure of 
the perturbative series [I] and essentially affects the definition of the renormahzed operator, as we shall 
see in a while. 

Now we want to consider the Liouville field theory as a perturbation of the free boson theory: 

5oM - j d^x ^^g^, <I>p„, = e^^ g = ^i. (3.5) 

Indeed, in the free boson theory the Virasoro algebras are given by the generators £° , £° obtained from 
the energy-momentum tensor components 

Joj = -2^t0^ = -\{d^vf- (3.6) 

With respect to this energy- momentum tensor the exponential operators Va{x) are primary as well, but 
they possess the conformal dimensions 

A^ - -a". (3.7) 
In particular, the conformal dimension of the perturbation field is 

Apert = A° - -b". (3.8) 

Consider the operator DmnDmnVa{x) in the vicinity of the point a = amn- Precisely at this point it 
vanishes. It is easy to show that in the framework of the free boson field theory 

DmnD,nnVa{x) = 0{{a ~ a,„„)^) (free boson). (3.9) 

Indeed, the free massless field ip{x) — ip{x^,x^) can be split (up to some zero mode, which does not 
interest us here) into the sum ip{x^,x^) — (p^{x^) + ip^(x^). Hence, the operator DmnDmnyaix) 
factorizes into the product (i'mne"'''"*-^ •')(5„j„e'^'''+'-^ •*). Each of the two factors is an exponential 
Qaip± niultiplied by a polynomial in the variable a. As far as it vanishes at a = amn it is of the order 
0{a — amn) as a — ^ amn, so that the product of two factors gives rise to 
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Taking into account (I3.9p . it is natural to define the operator 

^a\mn{x) = (3.10) 
yd CLran } 

which remains finite in the free boson theory in terms of derivatives of the field f at the point a — a„in ■ Its 
conformal dimension is equal to AjJ|^^ — —a} + mn. Its n folded fusion with the perturbation operator 
Vb generates operators from the Fock module over the exponential operator Va+nh of the conformal 
dimension A^^^^^ — (a + nfe)^. It is easy to see that the resonance condition 

A0|„,„ = A0+„, + n(l + (3.11) 

is just satisfied as a = Omn- In the vicinity of this point 

A-AJ,°^(A) = 2n5(a-a„„,) 

so that we may substitute residues in A by residues in a. We see, that the higher equation of motion 
(|2.17p takes the form of a resonance identity: 

Res $a|mn(a:) = 2B„„Kr,, _„ (x) . (3.12) 

a—ajnn 

This trivial fact will be the basis for obtaining resonance identities for the sinh- and sine-Gordon models 
below. 

The expansion (|3.4p gives rise to 

$a|™n(x) = + $„.„(A; X) + 0(a - a™„), (3.13) 

providing the definition of the renormalized operator A remarkable feature of the Liouville theory 

is that there is a nice way to get rid of the ambiguity. We can benefit from the transformation 

if-^ip + C, ^i^fic^''^. (3.14) 

Suppose the parameter A to be proportional to jj^ . We know that under the transformation p. 141) the 
1. h. s. transforms as 

^ a\mn ' ^ a\Taiv^ i 

while the first term in the r. h. s. as 

n ^4n&(a-a„„)TA , R \ 4nb(o-o„„) (a-45Mf)^(a-a„„))^ 
-^mn/\ ^ a+nb ' ^rnn^ ^a+nb^ 

Comparing both sides of (j3.13p we get the transformation rule 

$,„„(KAi^a:) ^ ($„,„(At/e-^^«;a;) +8e<5n62B,„„K„,-n(x))e'^'""«. 

The simple transformation rule $rnn — > ^mne"'""'' only restores in the case 6 = 0. Thus, A is dimension- 
less and may be set to 1. The operator 

^'mn(a;) = l'™„(l;a;) (3.15) 

with the transformation rule 

*mn(a;) -> *„„(a;)e"'""« 

will be considered as the uniquely defined renormalized form of $a„„|mn- 

Now, let us make sure that this definition is consistent with the common sense. For this purpose 
explicitly consider the simplest case m = n = 1. In the free boson field theory the field $a|ii is given by 

^a\ii = d-ipd+ipe"'^ (free boson), (3.16) 

since d-d+ip = 0. Hence, its value at the point a = an = is well defined: $o|ii = d-(pd+(p. Our aim 
is to define the corresponding operator in the Liouville field theory. 
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In the Liouville field theory the expression is not so simple: 

<^>a\ii=a-^d-d+c'"^ = {d-ipd+^ + a-^d-d+<fi)c'''^ = a-^d-d+ipe'"^ + d-ipd+'fi + 0{a). (3.17) 

According (j3.13p . p.lSp we have 

$„|n(a;) = a"^ • 27r^&e("+^)'^(^) + + 0{a). (3.18) 

It can be elicited from the perturbation theory that for small values of a 

a_5+^e'"^ = 27r/i6e(''+^)'^ + ©(a^). 

Hence, by comparing p.lSp with p.l7p we immediately obtain 

"^uix) = d^ipd+ifix). (3.19) 

This is the point, where the appearance of C instead of O'^, in the expansion (13. 4p is crucial. If we put 
there OJ,, the renormalized operator would have the form d-ipd^ip + ipd-d+ip. This is hardly could be 
considered as a renormalized form of the operator p. 161) at a = 0, since the second term is not invariant 
under the transformation (I3.14p and cannot be cancelled by adding e*"^. Besides, the operator could be 
extracted from the a-expansion of d+d-e"''^ and, hence, it would mean an unnatural reduction of the 
space of local operators after the perturbation. 

Now we want to rederive p.l9p in a more conventional way. For this purpose let us slightly move 
apart the points where the (p fields are placed in the product d-ip d+ip(x). Define the renormalized 
product by subtracting the pair correlation function; 

d^ipd+ip{x) = {d-(p{x')d+ip{x) - {d-ip{x)d+Lp{x)))\x'^x- 

Let us rewrite it as 

d-(pd+(p{x) = {d-ip{x')d+ip{x) + {ip{x')d^d+(p{x)))\x>^x 

= {d-ip{x')d+ip{x) + ip{x')d-d+ip{x))\x,^x - {'p{x')d-d+(p{x) ~ {'p{x')d-dj^'p{x)))\x< 

(3.20) 

The first parenthesis in the last line is the second order of the a expansion of d-d+e"''^ . It can be rewritten 
as 

{d^v{x')d+^{x) + ^{x')d-d+^{x)%,^x = a-^d-d+ {e-'^ - a(^)L^o 

= {a-^d^d+e^^ - a-' ■ 2nf,be'"^)\^^^ = (<f ,|n(x) - a-' ■ 27rt,be'"^^^%^^ . 

We used the equation of motion and the definition (j3.17D of $ci|ii. In the second parenthesis of p.20p 
we can apply the equation of motion: 

{ip{x')d^d+ip{x) - {ip{x')d^d+ip{x)))\,,^, = 2^/.6(^(x')e'"^(^) - (^(x')e^^(-)))|^,^^ 

= 27rM6^e^-=27r^6Ae(^+").|^^^. 

Here the second and third equalities are two reasonable definitions of the operator ipe'"^: as a limit of 
a separated product and as the a derivative of the exponent. We assume these two definitions to be 
equivalent. 

After subtracting two last expressions we obtain 

d-pd+ip{x) = ($,|n(x) - • 271^^66^^ - 27Tfib£e^''+''^'^)\^^^ 

= ($,|n(x) - a-' ■ 2nf,be(''+^^^<^-%^^ = ^nix). 

Though the manipulations are not quite rigorous and rely on some assumptions, the answer pays for it. 
We see that the construction in terms of resonance expansion conforms the common sense by providing 
reasonable regularizations of the operators defined by the perturbation from the free boson field theory. 

Finally, we succeeded to get rid of the ambiguities in renormalized operators in the Liouville field 
theory. More precisely, we swept the ambiguity under the carpet. It became possible due to the fact 
that the Liouville theory is a conformal field theory itself and we succeeded in non-dimensionalizing the 
constant A. We shall see below that in the case of massive field theories the analogs of the operators 
are intrinsically ambiguous. In particular, their vacuum expectation values contain logarithms of 
the dimensional coupling constant, so that a dimensional regularizer like the parameter A defined above 
becomes necessary. 
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4. Resonances in the sinh-Gordon model 



The sinh-Gordon theory is the field theory with the action 

S-shcM - J (^%^ - ch6^^ . (4.1) 

It can be considered either as a perturbation of the free boson field theory by the operator chbip or a 
perturbation of the Liouville field theory by the operator e"*"^. We shall need both points of view, but by 
default we shall use the former. Each time we use the second point of view, we shall declare it explicitly. 

Since there is no invariance like (j3.14p in the sinh-Gordon theory, the theory is non-conformal and the 
parameter /x is a genuine dimensional parameter of the dimensionality (mass)^"'"^^ . In fact, the theory 
is massive, so that the parameter M defined in (12. 7p is proportional to the mass m of the particle [5]: 

V^Q5^ (4.2) 

r(i^)r(H) 

Due to the loss of the conformal invariance the local operators can acquire nonzero vacuum expectation 
values (VEV). The VEVs of the exponential operators 

Ga = M-^'^'Ga = {Va{x)) = (vac| K (a;) | vac) , (4.3) 
are known exactly |9l[T0]. It is convenient to write down the factor Qa in the form 



g^^e-^'^-'^' cos^ 



X 



IT ^ N 7 N 7 { 7 0<p<l, 

aipfc\T^/l|aip/£\T^/l a, (l-p)fe \ ( 1 i a , (l-p)fc \ 



2'Q^2i^l2 Q ^ 2 ) \2 ^ Q ^ 2 



(4.4) 



with 7b being the Euler-Mascheroni constant. This form makes it apparent that the VEV as a function 
of a possesses no poles, while its zeros, being simple for generic values of 5, are located at the points 

Zeros: a = ±am„, m > n > 2, ne2Zorn>TO>2, m £ 2Z. (4-5) 

The VEV as an analytic function satisfies the equations |TOj : 

Ga — RaGq-a, Ga = G-a- (4.6) 

Below we shall need the notion and elementary properties of form factors of a local operator. Consider 
a model with one particle, like the sinh-Gordon model, on the plane. Let . . . ,0Ar) be the stationary 
state of N particles with the rapidities 9i > ■ ■ ■ > 9^ , s,o that |vac) — \0). Then the form factors of a 
local (or quasilocal) operator 0{x) are, by definition, the matrix elements 

(0;,..., 0^,10(0)101,..., 

In principle, for an integrable model the form factors can be found exactly, up to an unknown common 
factor, by solving a system of linear functional equations [S]. The words 'in principle' mean that there 
is no straightforward way to identify a solution to the equations with some particular local operator in 
the Lagrangian approach. 

For the exponential operators the form factors were found in [111112] . The common factor is fixed by 
the VEV, so that 



J-^,...,Wj^ 



Va{x)\9i, 9m) = GafaiOi, 9m, 9'j^ + ITT , . . . , 9[ + iir), (4.7) 



where fa are analytic functions of the rapidities and /a(0) = 1. To distinguish these functions from the 
full form factors we shall call them relative form factors. The relative form factors are dimensionless and, 
hence, /i- independent. It is important to note that in the case of the sinh-Gordon model the relative 
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form factors are finite for all values of a, so that the exponential operator Va just vanishes at zeros of 
the function Ga'- 

Vk+2,k+2+2l{x) = Vk+2+2l,2+2l{x) = V-k.-k-2l{x) = V-k-2l,-k{x) =0, fc, Z = 0, 1, 2, . . . (4.8) 

Besides, in contrast to the VEVs, the relative form factors depend on b meromorphically on the whole 
complex plane. 

Later we shall also need the notion of relative form factors extended to descendant operators. For 
any level [L, L) descendant (never mind, in the Heisenberg or Virasoro sense) V]^{x) of the exponential 
operator Va{x) the form factors can be written in the form 

{e[, . . .,0'j,\v^{x)\9u . . .,0m) = G^ff{0i, ■ ■ .,eM,o'j, + iTT, . . .,e[ + iTT). (4.9) 

Here is any normalization function proportional to (mass)^^°+'^+^. It makes the relative form factors 
dimensionless. We shall try to choose the form of the normalization factor in such a way that it took over 
as much nonanalyticities in b as possible. Below we shall make use of the conjecture that the relative 
form factors are analytic in b providing physically relevant answers for both real and imaginary values 
of b. One of the convenient choices of the normalization factor is M'^^^^^Ga, but we shall also consider 
other possibilities. 

The sinh-Gordon model can be considered as a perturbation of the Liouville model in two ways. 
First, we may consider the term with e^'^ as a part of the Liouville action, while considering the term 
with e~^'^ as a perturbation. Second, we may vice versa consider e~^'^ as a part of the Liouville theory 
and e*"^ as a perturbation. The two perturbation theories are related by the </? — — symmetry of the 
model. The Lp —ip invariance imposes some property on the form factors of the operators. For the 
exponential operators it looks like [TT| : 

{9[, . . .,e'^\v^a{x)\9,,. . .,0m) = {0[, ■ ..,e'r,\Vaix)\0i, . . . (4.10) 

The second of the relations (j4.6p is a particular case of this property. In other words, the particle 
creation/annihilation operators [13 are odd in the field tp. 

The first Liouville perturbation theory imposes the identification of the operators (j2.5p on the sinh- 
Gordon theory. On the other hand, the second perturbation theory imposes another identification, which 
is obtained from (|2.5p by the substitution a — >■ —a, Q — > — Q. Hence, the exponential operators of the 
sinh-Gordon model satisfy the identities 

Va{x) = RaVg^aix), (4.11a) 
Va{x) = R-aV.Q.aix). (4.11b) 

By using the properties of the VEVs (|4.6p we can rewrite these identities in the form 

Ga'Vaix) ^ GaVa'ix), if fl' - a G 2ZQ Or fl + a' G (2Z + 1)Q. (4.12) 

Note that the refiection property for the descendants (|2.1ip imposes relations similar to (14. lip on the 
descendant operators of the sinh-Gordon model as well. 

As it was shown in |14j the form factors obtained in [121 for the exponential operators are consistent 
with these identities: 

fa{0i,...,eN) = fa'i0i,.--,0N), if a' - a G 2ZQ Or a + fl' G (2Z+1)Q. (4.13) 

This symmetry property affects essentially the interpretation of the analytic function Go- In fact, it can 
only be interpreted as a vacuum expectation value on the strip 

- Q/2 < Rea < Q/2. (4.14) 

On this strip the exponential operators are normalized by the condition of the form 

{Va'{x)Va{x)) = Rq I \4aiQ-a) as |a; - a;|^ 0, a^asignRea. (4.15) 
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The first two terms in the numerator are sensible only on the boundary lines Re a = ±Q/2. The points 
a = ±Q/2 = Ta22 are exceptional, because Rq/2 — —1 and, hence, V±q/2{x) = 0. It is easy to check 
that the true primary operators at these points are ^±q/2(^)- other zeros of the VEV lie off the 
strip, where there is no reasonable normalization condition. We may think the off-strip operators to 
be just defined by the reflection equations (j4.1ip in terms of the on-strip operators. The zeros may be 
considered as a result of the unhappy normalization, though there is no way to improve this condition 
so that zeros would disappear without loss of the analyticity of the VEV. 

Another issue is related to the fact that the reflection function Ra possesses poles at the points 

Poles of i?a: a = afco, a — aok, fc = 1,2,3,... (4.16) 

The poles ako with k < 1 + and aofc with fc < 1 + lie on the strip < a < Q/2. This means that the 
short-distance correlation function {Va{x')Va{x)) for these values of a rises faster than \x' — a;|-4a(Q-a) 
as \x' — x\'^ — > 0. These fields can be normalized by the condition (see Appendix Ej) 

2(2 loff \x' — 

{Va{x')Va{x)) = for a = ±a,o or ±ao,, fc = 1, 2, 3, . . . (4.17) 

The corresponding fields in the Liouville theory are logarithmic. The existence of such fields was noticed, 
in particular, in |15| . 

Due to finiteness of the VEV Ga it is safe to say that there is no resonances for the exponential 
operators themselves. For even values of m it is easy to check that whenever the resonance condition 
for an exponential field is satisfied, the operator in the right hand side of the resonance equation is a 
descendant of one of the vanishing operators (14. 8p . The absence of resonances for odd values of m means 
that the corresponding perturbation contributions vanish. In the next section we shall see that this fact 
plays an important role in the sine-Gordon theory. 

Let us try to find the resonance of a level {N, N) descendant of the exponential operator Va (x) 
with a level {N',N') descendant of the exponential operator Va+(k-i)b{x) in the (fc + Z)th order of the 
perturbation theory, k,l = 0,1,2, . . k + I > 0. It means that Va+(k-i)bix) must appear in the operator 
product expansion of the initial descendant operator with k instances of the e^'^ operator and / instances 
of the e'^"^ operator. We have the resonance condition in the form 

Al + N^ A0+(,_;), +N' + {k + 0(1 + b'). (4.18) 

The solution to this equation reads 

lk + l + N' - N , 1 fk + l 



^ 2 k-i ^ +n^"'+r- ^'-''^ 

Note that this immediately discards the case k — I as a. resonance one. 

We are not planning to consider all this large selection of possible resonances, but specialize to some 
subclasses. First of all, let us restrict ourselves to the values a = amn with positive integers to, n as it 
was defined in (|2.14p . In this case k + I and N — N' must be divisible by /c — Z and wc obtain 

k + l = {l-k){l-k + n-l), (4.20a) 
N -N' = {l-k){l-k + n-m). (4.20b) 

Consider the resonances inherited from the Liouville theory. These resonances are related to the 
operators $a|mn as a — > For these operators N = mn, which admits a set of solutions to the 

equations (j4.20p . The solution k — I = n gives I = 0, N' — and corresponds to the contribution of the 
same field Vm.-n as in the Liouville theory. To find other solutions let 

I ^ k — m — s, 

which corresponds to the operator Vm.n+2m-2s- Then we can impose the conditions 

N' = s(m + n - s) > 0, 
, ('T^ - s){m + n-s-2) ^ {m - s){m + n - s) ^ .^21) 

2 ~ ' 2 ~ ' 

k,l£Z, k + l>0. 
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From these inequalities we immediately obtain 



< s < TO, TO — s G 2Z or m + n — s G 2Z. 

Besides, from (j4.5|) we see that Vm,n+2m-2s = for to G 2Z. Since N' < m{n + 2to — 2s) we conclude 
that the corresponding descendant vanishes as well and does not contribute the pole. Hence, nonzero 
contributions may only appear for odd to. 

Now recall that the contents of both the Liouville and the sinh-Gordon theories does not change, if 
we substitute 5 — > 1/& for a fixed value of the parameter M (or of the mass to) [TnilS]. This means that 
every nonzero contribution to the resonance residue must be reproduced in the perturbation theory after 
this substitution. It means that if we have a nonzero contribution from a level (N' , N') descendant of 
V„i^n+2m-2s, there must be a nonzero contribution from the same operator in the ch.b~^ip perturbation 
theory. But in the latter, the nonzero contributions can come from level {N',N') descendants operators 

Ki+L-2s',«' where 

N' = s'{m + n- s'), < s' < n, n - s' G 2Z or to + n - s' G 2Z. 

To reconcile both perturbation theories we have to require each operator V^''^„+2m-2s that give a nonzero 

contribution to coincide (up to a constant factor) to an operator V^_^2n-2s' «■ There are two possibilities. 
The first one is s + s' = to + n, since 

am,n+2m-2s ~ am+2n-2s' ,n ^ {n ~ s')b^^ + " "^)6 G 2ZQ , if TO + n = s + s' and s - TO G 2Z, 

but this is not the case due to its inconsistency with the conditions s < m, s' < n. 

The second possibility, which is consistent, is s' = s. It imposes some restrictions to the values of s. 
First, s must be less than both to and n. Second, 

am,n+2m~2s + a,„+2,i-2s,ri = (1 - TO - « + s)Q G (2Z + 1)Q, if TO + n - S G 2Z, 

(s) (s) 

and the operators n+2m-2s K^+2n-2s n '^^^ only coincide for even m + n — s. 

Besides, the operator n+2m-2s vanishes for even to, while lvi+2n-2s n vanishes for even n. Hence, 
these operators can only make a nonzero contribution when both to and n are odd. 

Let us summarize the results. There are two types of resonance contributions: 

■ *a„„.|mn Vm.-n, k = n, 1 = 0] 

: *Q™„|m« yrn!n+2rn-2s ^ith N' , fc, I from (H?^ and < s < TO, 71, .s G 2Z, if TO, n G 2Z + 1. 
Here the sign ^ means 'is in resonance with' and v}^f,^^2m-2s some level {N',N') descendant of the 

operator Vm,n+2rn-2s- 

The R_(- resonance is the same as in the Liouville theory and, since it originates from just one term 
in the perturbation expansion and this term is the same as in the Liouville theory, it is perturbatively 
exact and the coefficient at the field Vm,-n coincides with that obtained in [3]. The Ri resonances are 
specific for the sinh-Gordon model, but they are perturbatively exact as well. 

To write all the resonance identities in a compact form it is convenient to define the sets 



{s G 2Z I < s < TO, n} for to, n G 2Z + 1, 
otherwise. 



Then the higher equations of motion in the form p. 121) are modified as follows 

Res <i>,|,„„(x)=2i?„„K„,-„(x)- '2C^lyi%2rn-2si^), (4-22) 

n.—n. ' ' ' ' 



sGZ„ 



where Cmli — A*'™ * ^^Cmli are some coefficients. The numbers Cmli depend on the normalization 

(s) 

of the operators n+2m-2s- Since the explicit form of these operators is unknown, they are arbitrary 
except for the case s — Q, where we may identify V^\j^2mi^) — Vm,n+2m{x)- 
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For even m or n the resonance equation retains the form of the higher equations of motion of the 
LiouviUe theory: 

Res <^>a\,nn{x) = 2B,nnV,n,^n{x) , if TO G 2Z Or n E 2Z. (4.23) 

a—amn 

For odd m and n the cases m — 1 and n ~ 1 are the simplest ones. We have 

Res ^a\mn{x) = 2BmnVm-„{x) ~ 2CmnVm,n+2mix) , if TO, Tl G 2Z + 1 and TO = 1 Or n = 1, (4.24) 
a—amn 

where the coefficients Cmn = Cmn can be found exphcitly. In principle, we could find them by taking 
appropriate multiple integrals coming from two-point correlation functions, but it is possible to circum- 
vent explicit integrating. We can make use of the fact that in this case the descendant level mn of 
the operator $a|,„„ is odd. On the other hand, for a generic value of a the Fock module coincides 
with the Verma module of the Virasoro algebra. Accordding to the widely accepted conjecture (see 
e.g. [13), it is spanned by the vectors generated by the commutative integrals of motion I_2A;+i of odd 
spins 2fc — 1 and by the Virasoro elements £_2fc of even spins 2k. The action of the operator I_2fc+i 
on an operator 0{x) is just the commutator with the integral of motion /_2fc+i in the space of states: 
X-^k+iO^x) = [0{x)^ /_2fe+i]. Since the vacuum is annihilated by all commuting integrals of motion, we 
have {X-2k+iO{x)) = {[0{x), I^2k+i]} = 0. As far as every Virasoro descendant of an odd level can be 
represented as a linear combination of commutators of integrals of motion with other operators, its VEV 
is zero. Hence, ( Res a=a„„ ^a\m7i{x)) = 0, whence it follows that the VEV of the r. h. s. of (|4.22p must 
be equal to zero. The last is satisfied, if 

Gm. n B.,nn J^"'' " ; if TO, n G 2Z + 1 and to = 1 or tt. = 1. (4-25) 



Since am,n+2m + 0'm,-n = (1 ^ m)Q E 2'LQ we obtain from (|4.12p that 

Vm,n+2ni[X) = — if TO G 2Z + 1. 

Hence, we can rewrite the equation (|4.24p in the form 

Res ^almnix) = 2B„„(K„ _„ (x) - V-a,„ ^„{x)), m,n G2Z + 1, TO 1 or n = 1, (4.26) 

a—amn 

which nicely generalizes the standard equation of motion of the sinh-Gordon model corresponding to 
TO = n = 1. Some checks of (I4.23|) . (j4.26l) within the form factor approach are given in Appendix IB] 

Now we are ready to define the operators "^mnix) in analogy with those in the Liouville field theory. 
It is convenient to write the dimensional parameter A, defined in p.4p . as kM, where k is an arbitrary 
dimensionless constant. Then we get the expansion 

<i>a|rnn(a:) = — ^ ( i?„,„ (k+M)4"''(«--'"") K+„b (x) 



^ c(:)(^L^)m)-4(™-)''(-'^™.)t/(^)_^^^(^)) 



sez„ 



^'™„(k;x) + 0{a - a 

mn 

), (4.27) 



where k denotes a vector made of all k variables; here k = (k_|_, k!^'', kI^-*, . . . , k?^'"*-™"^'"""'^''-'). Let 
1 = (1,1,...,!). We shall denote 4'mn(l; x) = 5'm„(a;). Now we want to rewrite this expansion in terms 
of the relative form factors. Define the relative form factors /*„j„, fmn of operators ^a\mn, ^mn by 
the equations 

(vac|$,|„„(O)|0l = A"+'"<'^~'"""^/^G,+„,/*_(0l, 

J (4.28) 
(vac|vI'™„(«:;O)|0) = A"G,n,-„/*i^H^), 



where is a short notation for the sequence 6i, . . . ,9pf of rapidities of N particles. Besides, define the 
relative form factors of the operators v''-''\, . as 

(vac|eL-.).(0)l^^> = M^^("+"-)G,_(,„_.),/fl) ,(^1. (4.29) 
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The prefactors are chosen in such a way that the relative form factors were dimcnsionless and finite. 
By substituting them into (|4.27p we obtain 



+ /,L(^) + 0(a-a™„)- (4.30) 

For even m or n the expression is simple: 

/aVW = ^^"1^°^"^^^^ + fZniO) + 0{a - a_), (4.31) 
while for odd to, ri and to = 1 or n = 1 it is not much more complicated: 

-f* fQ\ 2S,„„(/a+„fc(0) — Ga+nb^a-mbfa—mb{d)) 5, , , , , 

fa\mnW = — " + fmrM + 0[a - a„„). (4.32) 

The K;- dependence is given by the rule 

- m^~s)CitlS;U5a-(m-s)bU.n+2n^~2siO)logK^ll (4.33) 

In the operator language the transformation rule reads 

*m«(K;a;) = *mn(2^) - 8nbB„inVm, ~n{x) log K+ 

- J2 mm~s)Ci^lVm.2m+n^2s{x))\0gK.'-_!\ (4.34) 

This transformation rule incorporates the logarithmic ambiguity of the definition of the renormalized 
operator. From the form factor point of view the operator \l/(a;) looks most natural, but it is necessary 
to stress that all operators "^{k'jX) are, in principle, equitable. Some particular values of k can be 
distinguished by some particularly nice properties. For example, if to or rt is even, the choice log k+ = 
/*„(0)/8n6;B,„„ makes the VEV of the renormalized operator vanish. 

5. Resonances in the sine-Gordon model 

Formally, the sine-Gordon model with the action 

Ssg[^] = J d^x (^^^^ + 2(-m) cos/3(^^ , < /32 < 1, ^ < 0, (5.1) 

can be obtained from that of the sinh-Gordon theory by the analytic continuation to imaginary b = — i/3 
for fixed fl. In terms of the parameter p this theory corresponds to 

p<0. 

Surely, we want the mass to to be a real positive number. Then according to (|4.2I) the parameter M 
must be a complex number with the argument m{p — 1). This is consistent but inconvenient. Thus we 
make the substitution 

so that the new parameter M be real. We shall mean by M this real number in the rest of this section. 

Some quantities can be indeed obtained by the analytic continuation. Nevertheless, there are some 
essential differences between the two models. Physically, these differences are the result of different 
particle spectra. The sinh-Gordon model contains the only particle, which can be identified with the 
boson ip. The sine-Gordon model for any admissible values of (3 contains a pair of topological solitons 
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— kink and antikink, which can be considered as elementary particles. For < 1/2 they form a 
series of bound states — breathers, the lightest of which, the first breather, can be identified with the 
'continuation' of the sinh-Gordon particle. 

We retain the notation Va{x) for the exponential operators, though for physically admissible operators 
the parameter a is purely imaginary, a — ia. The degenerate points amn correspond to 

- l—H^R-^ - 
Otmn 2 2 

The form factors of the exponential operators can be again factorized into the VEV and the relative 
form factors: 

(vac|K(O)|0l£l,...,0JVeAr) =Ga/a(ei,...,0Jvk...e„, /a(0) = l, (5.2) 

where the variables £i = +,—,1,2,..., — 1] denote the sorts of particles, + for the kink, — for 

the antikink, 1,2,... for the breathers in order of increasing mass. The relative form factors that only 
contain the first breather are obtained from those of the sinh-Gordon model by the analytic continuation 
in the parameter b: 

/a(^l,...,^jv)l...l= /a(ei,...,f?)|,_i^. (5.3) 

Moreover, we conjecture that the first breather relative form factors for any local operator can be 
analytically continued from the sinh-Gordon case. 

Nevertheless, since the kink form factors do not satisfy any reflection conditions like (j4.13p . the 
equations (j4.1ip do not hold and the operators Va, Vq^a and V-q-a are essentially different. 

The vacuum expectation value Ga = {Va) is similarly represented as A/^^° Qa^ but the factor Qa for 
the sine-Gordon model cannot be obtained by analytic continuation of that in the sinh-Gordon theory. 
The obstacle is the fact that the set of zeros (|4.5p become everywhere dense for negative p in generic 
position. It was conjectured in [TU] that the VEV has the same form in both cases being written in the 
integral form, but the analytic structure of this integral in both regions is completely different. Rendered 
to the form of an infinite product of gamma functions the factor Qa reads as 

g ^e-^^^'' TTe^ ^ ^ ^ - ' ^ 5^ ^ «<0 (5 4) 

ii rMi-f)r(i-§ + ii^)r(i + § + ii^V ^ ' 



This function also satisfies the equations (|4.6|) . but already possesses poles. It is easy to list zeros and 
poles of this function; 

Zeros: a — ±amm m > n > 2, mS2Z-|-l, n G 2Z; 

(5.5) 

Poles: a = zLamn, m > n > 1, me 2Z, n G 2Z + 1 or m > 2, n < 0, m G 2Z, n G Z. 

The interpretation of the zeros is quite different from the case of the sinh-Gordon model. In that case 
the relative form factors fa{- ■ ■ ) were finite, and a zero of the VEV meant the whole operator vanish. It 
was the result of an unhappy normalization of the corresponding exponential operator off the 'physical' 
strip (I4.14[) . In the sine-Gordon model the normalization condition 

{V^a{x')Va{x)} ^ I , ^ 14^0 as|x'-x|2^0 (5.6) 
\x' — xp"^" 

is not so vulnerable as ([4.15^ : it makes sense for any value of a. It means that no exponential operator 
can vanish, and, hence, if the VEV of the operator Va for some particular value of a vanishes, some of its 
relative form factors tend to infinity at that value, so that the corresponding absolute form factors remain 
nonzero. The relative form factors that only contain breathers are always finite. Hence, the infinities 
must appear for the soliton (kink-antikink) relative form factors. Lukyanov's integral representation |12) 
for the soliton relative form factors of the exponential operators is convergent on the strip 

\lma\ < - P/2 (5.7) 

and diverges as a — > ±i(/?^^ — 13/2) = ±032. The values a — ±032, which bound the strip, are, indeed, 
the closest to the real axis zeros of the VEV. Every other zero of the VEV lies off the strip (|5.7|) all the 
more. This explains how an exponential operator can occur nonzero even if its VEV vanishes. 
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The existence of poles suggests the presence of resonances for exponential operators. To verify it 
consider the conditions (|4.20l) for = 0. Let I — k~m — n — s. Then k + I = (m — s — l)(m — n — s) 
and, hence, 

(m ~ s ~ 2)(m — n — s) , (m — s)(m — n — s) , . , , 

k = - -, l = - N' = s(m-n - s). 

It can provide resonances for m > n, m > 2, if either m — s G 2Z or m — n — s G 2Z and < s < m, m — nl^ 
This is possible in both cases of even and odd to. But for odd m in the sinh-Gordon model there are 
no poles even in spite of satisfied resonance equations. Hence, we are sure that the corresponding 
perturbation integrals vanish, and there is no poles in the sine-Gordon case as well. Hence, we limit 
ourselves to the case of even to. We get the resonances 

Vmn ^ l^m,2m-n-2s' 0<S<m — 1,TO — 71, SGZ, ifTO>2, n<TO, TOG 2Z, 71 G 2Z + 1, Or 

< s < m — l,m — n, sG 2Z, if to > 2, 77 < m, to, 71 G 2Z. (5-8) 

Here again VF^''_2„_„_2s is an (iV', N') descendant of Vm.2'm~n-2s and W^^jm-n = Ka,2m-n- We expect 
resonances not only for odd positive 77, as it can be seen from (I5.5p . but for even positive 77 as well. This 
suggests that the finite values of VEVs for even n do not mean that all form factors are finite: the kink 
form factors must possess poles. In other words, we may describe the poles and zeros of Ga as follows: 

'Physical' zeros: a = ±amn, m > n > 2, to G Z, 77 G 2Z; 

(5.9) 

'Physical' poles: a = ztomn, m > n, to > 2, 777 G 2Z, 77 G Z. 

We assume that at the 'physical' zeros the breather-to-soliton form factor ratios vanish, while at the 
'physical' poles the correlation functions of the exponential operator possess a pole due to the resonance 
phenomenon. The VEV of the exponential operator Ga is regular at every point that is a 'physical' zero 
and a 'physical' pole at the same time. 
The resonance identity reads 

Res Vaix)= J2 ^^lW^™,L-«-2., (5-10) 

where Dmli = /i^™~*~^'''™~"^*^I?mi are some coefficients and 

{{s G Z I < s < TO - 1, TO - n}, if TO G 2Z, 77 G 2Z + 1, TO > 2, 77 < TO, 
{s G 2Z I < s < TO - 1, TO - 77}, if TO, 77 G 2Z, TO > 2, 71 < m, (5-11) 
otherwise. 

Due to the relation a^n + o,m,2m-n = ijn + \)Q G (2Z + \)Q the breather relative form factors of the 
operators Vmn and Vm,2m-n coincide for even to: 

fmn{())g — fm.,2m-n(0)g, Ei = 1,2, .. . (5-12) 

This identity makes it possible to establish the coefficient at the operator Vm,2m-n- In the vicinity of 
the point a„in the VEV admits the expansion 

Ga = m2('^^"^"-^) lAl^ + G(°) + 0{a amn)] , GL"1 = M^^'^^g^^l (5.13) 

\ 0"mn J 

From (|5.12p and f|5.13p we get the proportionality coefficient 

D^l - J^^- (5.14) 



■^There is also a set of solutions for m < 0, but it can be obtained from this set by means of the substitution m — ^ 2 ~ m, 
n — > 2 — n, which corresponds to amn — > —a,mn- 
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Besides, the identity (|5.12p implies that the total contribution of the operators W^2m-n~2s '^i^h 
s ^ into the breather form factor is zero. Of com'se, their contribution to soliton form factors can be 
nonzero, so that the identity (I5.12p is unlikely to hold for solitons. 

To define the regular part consider the expansion 

Va{x)^ z ^-^^^^^^^^^^^ + V„U!^;x) + 0{a-a„,n). (5.15) 

Here v_ = {v'"'^^) seY„^„ are free parameters. 

To establish the form factor relations let us define the relative form factors g'^^^ of the operators VF^*) 

as 

{-^^\W^\^^,,-s)Mm = M'^^"^'-'^^'^^^^^ (5.16) 
Substituting it into (jS.lOp . we obtain 

Ql-n^Imn0),^gm,2m-n ^ 'D[:^n9t]2m-n-2s0)e- (5-17) 

Then from the regular part of (|5.15p we obtain 
(vac|l4„(z.;0)|^) = Gt]jmn{e)e + G[-^^ fl^^iO), 

seY„^,i 

+ {G'n.2m-n - 4^m,2m-n(m - 71 - s)6 log ^z^'') ) gt]2m-n-2si(^) t) , (5-18) 

where primes mean derivatives in a: • ■ ) = lafi^i' ' ' )|a— a ^^''^ "^^^ particular choice 

log iz^^-* — log vi"^ = 



4(n - TO - s)b ytJiTi, 2777-77 g^~,P J 

(s) 

is especially nice, since the first term is cancelled by the terms proportional to 9^ 2m-n-2s- 

(vac I K77„ (Z^; 0)1^)= Ct]^ fU„ {ff)e - Af-^-^™" ^777,2777-77 ^ ^™nffm!27n-77-2. 

seY^„ 

Below we shall need another choice: li = I- In the breather sector, due to (|5.12p . (|5.17p we have for 

^^77777 (■^) ^7777 (1_5 • 

(vac|V;„„(0)|^) = G(^if„Ue)r+ 2GL^)/;„(0lr, = G^^i - Gi-^^g;^],^_S,^^,^^„. (5.19) 

Now consider the mn level resonances. There are two cases: the case of odd to and that of even m. 
For odd to the VEV of the operator Vm,-n is finite, and the resonances are very similar to those in the 
sinh- Gordon case. It is easy to see that the only difference from (|4.27p is that the sum in the parenthesis 
must be taken over the set 

^, ^UseZ\0<s<m}, ifTO,nG2Z+l, 

[{.s e 2Z+ 1|1 < s < m}, if TO e 2Z+ 1, n e 2Z. 

The case of even to differs greatly. On one hand, the operators Vm,n+2in-2s do not vanish in this 
case, so that we have to admit an expansion like (j4.27p . but with summation over the set 

^, _ fls e Z|0 < s < to}, if me 2Z, ne 2Z + 1, 

~ [{s e 2Z|0 < s < to}, ifm,ne2Z. ' 
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On the other hand, the first term in (|4.27p contains the operator Va+nb, which itself has a pole at a = a„m 
according to (|5.10|) . (|5.11|) . since amn + nb = am,-n- Besides, this expansion must be consistent with the 
form factor identity (j4.3ip in the breather sector. 

First of all, notice that the ratios Ga/Ga+nb are meromorphic functions of b on the complex plane, 
whose explicit form can be extracted from (j4.4p for real values of b and from (|5.4I) for imaginary ones. 
Besides, for generic values of b they do not possess poles at a = amn for any m > 0, if n > 0. Hence, the 
numbers 

m,n>0, (5.22) 



are well defined: 

Kmn = (A/fe2(i-p))2(Qb-™)" 

^ TT ~ 26am^„+i_2fc)r(l + Q -'-a„i,Ti+l-2fc)r(^ — Q~"'"am,n+l-2fc) ^ 

k=\ ^ 26a„,„+i_2fe)r(l - (g-iarn,,„+i_2fc)r(i + (5~^am,„+i_2fe) 

Evidently, 

Kmn = 0, if TO > n > 2, n e 2Z orn > TO > 2, to G 2Z. (5.24) 

As a consequence, the definition of $a|„m in terms of the form factors is reliable in the sine-Gordon 
model as well. Thus we may continue the identity (|4.3ip analytically to the imaginary values of b and 
multiply it by the quantity /t"+^"'("~"™")/'3G'a+n6 for the sine-Gordon model, which possesses a pole at 
a = Omn- We obtain in the breather sector (e^ = 1,2,...) 

(vac|$,|_(0)|^) = /i"+2„(a-a„„)/Q / 2g„,„G,+„,/,+„fc(g),^ ^ G,+„,/,*„ (^1,-) + 0(1) 



0(1) (5.25) 



as a ^ amn- By comparing with (|5.19p . after generalization to the soliton sector, we may conjecture 



^mn) s^y 

, 2B™„ym _„(i/; a;) + *mn(a^) 



0(1). (5.26) 

(X Ojmn 

The finite operator ^ defined by this expansion possesses in the breather sector the form factors 

(vac|*„„,|0l) = A"GL:Ll/*„(elr. (5.27) 

Notice that the operator appears as a finite part of the expansion of (a — amn)^a\mn rather than 
of ^a\mn itself. It mcaus an important circumstance that the operator ^a\mn{x) needs a multiplicative 
renormalization beside an additive one. 



6. Conclusion 



We have used some consistency conditions to establish the form of the resonance identities and to define 
renormalized operators in the Liouville, sinh- and sine-Gordon theories. We found that in the case of 
the sinh-Gordon model the resonance identities for some operators remain the same as in the Liouville 
model, while for the other operators the resonance identities are modified getting a finite number of extra 
terms. These extra terms are perturbatively exact in the sense that for any given correlation function 
each of them appears in a single term in the perturbation expansion. We have seen that the sine-Gordon 
model possesses a much more complicated structure of resonances. The resonance identities are expected 
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to be exact operator identities that may help us to estabhsh the correspondence between the bootstrap 
form factors and local operators of the Lagrangian field theory. 

Strictly speaking, most of the expressions given here are conjectures. Their consistency to the form 
factor identities in some simple cases encourages, but we need a more rigorous derivation. An explicit 
computation of integrals in the perturbation theory would provide a rigorous check for these conjectures. 
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Appendix A. Proof of the normalization condition (|4.17p 



Here we derive the normalization condition at the poles of the reflection function Ra in the framework 
of the Liouville field theory, where we do not have to restrict the equations to the strip (j4.14l) . 

Let us start with a general remark. The pair correlation functions {Va{x)VQ-a{0)} , {Vaix)Va{0)) 
correspond to the points where the structure constants possess a pole [B]. It means that the two-point 
correlation functions in the Liouville theory formally diverge. Nevertheless, there is a way to define finite 
two-point functions as residues of some three-point functions. First define the functions 

(F5(oo)K(x)^Q-a-5(0)) = lim Ro^\xoo\^^' Res {Vs^{x^)Va{x)VQ^a-sm, 



.+s 



(A.l) 



(A.2) 



{Vs{cx)Va{x)Va+sm = Ra+5{Vs{(X^)Vaix)VQ.a-Sm. 

Then the two-point functions are obtained in the limit (5 — >■ 0: 

{Vaix)VQ^am = {V5{^)Va{x)VQ-a-5m\s^0, 
{Vaix)Vam - {Vs{<X,)Va{x)Va+5{0)))\s^O- 

It explains why we introduced the factor Rq^ into the definition of the residue correlation functions: it 
is necessary to conform the natural identity (Vb(x)Vb(O)) — 1. 

Now turn to the special values of a we are interested in. Let a, be any of the points ako or ooz, which 
are poles of the refiection function. Since Rq-a = it h^-s a simple zero at a,: 

RQ-a = -(a - «*)-R(3-a. + 0((a - a,)^) as a a,, 
where R'^ is the derivative of Ra- Now consider the product 



^b-a.-5^5MK.+5(0)|,^„ = - -RQ-a,-sV5{^)Va.+m 



6^0 



Since all the operators Va, (x) are finite (because they possess finite three-point functions with nearly 
all exponential fields), the operators Vg-a, (x) vanish. Hence, we may substitute the limit in the last 
expression by a derivative: 

^b-a. ^0(00)^^.(0) = - ±Vs{^)VQ^a,^siO) 



5=0 



By using the definition (jA.2l) we obtain 



-R'r 



dS 



{Vs{<x)Va,{x)VQ-a,-5ix)) 



5=0 



Taking the explicit form of the three-point function in the r. h. s. from (jA.l|) we see that 



{VaAx)VaM) 



-R'r 



'-a.-^-O 



4|^|2A.-2A„.-2A„.+. 



d6 



2a* 



5=0 



^Q-a.^0|xP-(Q-a.) 



log|a;|" 



which leads to the short range asymptotics (|4.17p in the sinh-Gordon theory. 
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Appendix B. Form factors of the level < 3 descendants and resonance identities 

In [14j a construction for breather form factors in the sine-Gordon model was proposedlf] First, we briefly 
review it in the Liouville/sinh- Gordon notation used in the present paper. 

Let ^ be a commutative algebra with the generators c_„, n = 1,2,..., and A be its copy with 
the generators c_„. Let be the associative algebra generated by elements of both algebras with the 
additional commutative relation 

, m ^ fo, m€2Z+l, 

[Cm, C-„J - 2 ^.^2 ^ ^rnn X ^ ^ 

To any element g of the algebra A"^ we associate a sequence of functions {P^{xi, . . . ,Xk\yi, ■ ■ ■ , yO}fe°/=o 
defined according to the following rule: 

P\X\Y) - 1; P^-{X\Y) = S,,{X) + (-)"-i5„(y); P''-{X\Y) = 5_„(y) + (-)"-i5_„(X); 
pki9i+k2g2^X\Y) ^ kiP3^{X\Y) + hP'^iXlY), h,k2eC, 91,92 eA""; 
P'''''{X\Y) ^ P''{X\Y)P'''{X\Y), h,h'eA; p'^'^' {X\Y) ^ p'^{X\Y)P''' {X\Y), h,h' e A; 
P'^''\X\Y) = P'\X\Y)P'^'{X\Y), he A, h' eA, 

where 

k 

X={xi,...,Xk), Y ^{yu...,yi), 5„(X) = ^x^ 

i=l 

It is important that in the last line the element h' precedes the element h in the product h'h. Let 
9 = h(E)h' e A(E) A. Define the element [g] = hh! e A? . The linear map [•] puts a 'barred' element to the 
right of the 'unbarred'. The form factor construction uses the functions where we first need 

to push the 'barred' elements through the 'unbarred' ones to the left and then apply the above rules. 
Define the constant A' and the functions R(6) and /(z) according to 

1 \ /■'^^ t 

„ ■ exp / dt^- 

2 sm ) Jo sm t 

°°dt shf sh^sh^ii^ \ 




(z + cj)(z-a; 1) _ . p 



z^- 1 

Then the relative form factors of some descendant V^{x) [g € A® A) oi the operator Va{x) are given by 

N 



a 1 

where j^{X) are the symmetric functions given by 



Here the sum is taken over all exact partitions of the (multi)set X into pairs of nonoverlapping sets X^ 
and X+. In the simplest case g = \ ® \ the functions fa{S) = fa^^{d) are the relative form factors of 
the exponential operators Va{x). 

Now we are ready to describe several particular cases, related to the resonance operators described 
in the paper. We can introduce some particular elements h[j\ Z = 1, 2, 3 on each level I of the form 

hi'^ = (B.l) 

COSTTZ^ 



construction for soliton form factors was also proposed, but it needs some refinement due to problems with the 
convergence of integrals. 
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" sin 7rp + sin ' " 

^(3) ^ Ji 1 3'=°^'^'-^ B.3 

cos 2ttp + cos 2ttv 

We use them to define the operators 

= gi')=ft«e. (B.4) 

It can be easily checked that the operator 



\Trm J 



possesses a simple pole at a = and 



ResT/W=-3)%+a_^, 

a=0 \TTm J 

where m is the mass of the first breather (j4.2p . The corresponding breather form factors of this operator 
indeed satisfy the identities [18: corresponding to the equation of motion 

d+d-tf — Air^bshbip = 4:TT{—fi)(3 ain/Sip, 

which is the resonance identity for ^a\ii from (|4.24p . It is easy to obtain the answer for the renormalized 
operator ^'11(2:) = d+(pd-(f{x) in terms of relative form factors of exponential operators: 

(vac|9+^9_^(0)|^) = i (E™-.«'*) (E"*-'°"'') f^{^)e^2l,^JibK^,\ae)e- fU0)e) 

- 2nfxbK{^^ ({L + mog^,+)Me)r+{L + 4MogK_)/_6(^)r) . (B.5) 



Here is the mass of the eth particle (so that nii — m), and L — G[i/Gb- It seems to be reasonable to 
choose K+ = K_ . For this choice the form factors for the states consisting of the first breathers vanish for 
odd number of particles. Intuitively, this corresponds to the fact that the operator is even with respect 
to </J. 

The most convenient choice is logK+ = logK_ = — i/46, since the expression in the last line com- 
pletely vanishes. For this choice the zero- and one-breather form factors vanish, while the two-breather 
one reads explicitly 

(vac|9+(^5_^(O)|0il,02l> = TT''p{l-p)\'m'R{ei - 62) ((e^o,-e2)/2 ^ ^(e,-e,)/2-^2 „ 2cos7rp) . 

In the limit b ^ 0, m ~ const it tends to its free field value 7rm^(e^i^^2 + c^^^^^). 

The following results have been checked up to the 4-particle form factor analytically and up to the 
8-particle form factor numerically. 

For Va (x) we have the conjectural identities 

Res V^^\x) = Va^ mn = 12,21, 

27r sin irp 



a—a^r 



which can be compared with the $a|i2 ^-nd $a|2i resonance identities (j4.24l) . This comparison gives 

V!l^\x) = j^^-"" ^a\mn{x) + 0{{a - a„„)"), mn = 12, 21. 
47ri?,„„sm np 

(2) 

Of course, we can extract the finite parts of the operator Va , but we cannot be sure that they co- 
incide with the corresponding operators ^12, ^21- They may contain contributions from C^j^C^j^Vmn, 

■C?-l^-2Knn, •C-2^?_iKnn- 

For V^-^'ix) we can also conjecture some identities. First, is was checked that the relative form factors 
only possess simple poles at a — 013, as 1. Then we have 

Res V}'Hx)= .^3^"" (K„._„(a;)-^_,,„_„(x)), mn = 13,31. 
dZTT sm np cos^ 7rp 



a—a„ 
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